function Fig3_b()
lb = 25.66 / 6^0.5;
EC = 56.1 * 6^0.5;

% This data is prepared by matlab code dataPrep_impurityA_multi_filling.m
% dataPrep_impurityA_multi_filling.m is located in .../processed data/
data_r = readmatrix('processed data/multi_filling_2D_r.csv');
nu = data_r(2:end,3);
n_nu = length(nu);
r = data_r(1,4:end)*lb;
phi_r = data_r(2:end,4:end)'*EC;



Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

%% make uneven plot
n_nu_u = 121;
nu_u = linspace(-1,0,n_nu_u);
for i = 1:n_nu_u
    [~,id] = min(abs(nu-nu_u(i)));
    phi_r_u(:,i) = phi_r(:,id);
end



%% plot

figure()
imagesc(nu_u,r,phi_r_u);

ax=gca;
ax.Color=Bcolor;
ax.XColor=Fcolor;ax.YColor=Fcolor;
xl=xlabel('filling factor');
yl=ylabel('r (nm)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;xl.FontSize=Fsize;
ax.YDir='normal';ax.XDir='normal';

cMap=getColor('Magma');colormap(ax,cMap);
ax.CLim=[-0.8,0.4];
xlimit = [-0.85,-0.15];
axis([xlimit(1),xlimit(2),0,max(r)])

end
